Rashba precession in quantum wires 
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The length over which electron spins reverse direction due to the Rashba effect when injected 
with an initial polarization along the axes of a quantum wire is investigated theoretically. A soft 
wall confinement of the wire renormalizes the spin-orbit parameter (and the effective mass) stronger 
than hard walls. Electron-electron interactions enhance the Rashba precession while evidence is 
found that the coupling between transport channels may suppress it. 
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Spin transport has regained considerable interest in re- 
cent years, partly because of the prospective for quali- 
tatively novel electronic devices and partly since funda- 
mentally new mesoscopic quantum coherence phenomena 
could be investigated in the magnetic degree of electronic 
freedom. A prime example is the spin transistor proposed 
by Datta and Das [1]. 

Here we investigate how Rashba spin precession [2] is 
affected by many body effects [3] as expected to be im- 
portant in semiconducting devices. We consider electrons 
moving along a one-dimensional (ID) quantum 'wave 
guide'. Also the cases of two interacting channels and dif- 
ferent transversal confinement potentials are addressed. 
Of prime interest will be the length Ar over which spins 
initially polarized along the channel axes reverse spin- 
direction while moving along the channel. 

In a 2D layer spin splitting is a consequence of spin- 
orbit coupling 



H so = a(a x p z - a z p x ) 



(1) 



proportional to the Rashba parameter a and thus to the 
intrinsic or by means of gates externally applied electric 
field in y-direction, perpendicular to the layer [2]. H so 
is also proportional to the momentum of spins (j> x ,p z ) in 
the 'active' x — z-plane, a x z are Pauli matrices. Preces- 
sion then occurs on the length scale |fc + — fc_| _1 of the 
spin split wave numbers k± at the Fermi energy. If the 
effective mass approximation e(k) = k 2 /2m is applicable, 
as for many semiconductors, the Rashba length 



27r/|fc + — k—\ = it/ma 



(2) 



does not depend on carrier density or Fermi energy. 

With electrons being confined in z-direction, to zero'th 
order in _ff so this picture carries over to ID. Only to 
higher order H s ° will weakly mix different transport 
channels of the wave guide which up to 0(a 5 ) can be 
accounted for by renormalizing a — > a* in e±(k) and, 



within the effective mass approximation, by m — > m*. 
In the latter case the Rashba length modifies according 
to (ma)" 1 — > (m*a*) . For a quantitative estimate of 
both renormalizations the intra-subband cigenfunctions 

i>kns(x,z) — e (f> n (z)(cos(maz)\s) +i sin(maz)|— s)) 

(3) 

are needed, which are plane waves of momentum k along 
the wave guide and, without inter-subband scattering, 
slightly modified subband states cj> n (subband index 
n) in z-direction, spin polarized s — ± on the axes at 
z = 0. For a harmonic confinement (subband energy 
ujq), as relevant for example for samples defined by gat- 
ing [4], the perturbative estimate yields a* — a(l — 77) 
and m* = m(l + 8r/ 2 ) in the ground subband. The 
dimensionless parameter 77 = (wma/2) 2 measures the 
bare value of the spin-orbit strength a in (1) by com- 
paring the wave guide width w = 2/^/mwo with the 
length (ma) -1 . For a hard wall confinement on the other 
hand (again of width w), as possibly more relevant for 
wires fabricated by the cleaved edge technique [5], the 
renormalizations become a* — a(l — (1/6 — l/ir 2 )r]) 
and m* = m(l + 3(4/37r) 6 ?7 2 ) , i.e. they are signif- 
icantly reduced compared to the soft wall case since 
1/6 - 1/tt 2 « 0.065 and 3(4/3tt) 6 /8 w 0.002. 

The main idea of the transistor operation [1] is to vary 
the strength of the electric field and thereby it /ma. This 
is achieved by a charging a gate parallel to the layer. 
Unless carefully compensated by a second back gate [6], 
however, the carrier density in the structure will change 
at the same time. Ignoring interactions this would be 
unimportant when e(k) = k 2 /2m* but the interaction 
strength depends, even quite sensitively, on carrier den- 
sity through the r s -parameter. In order to investigate 
possible consequences for the Rashba precession we em- 
ploy the Tomonaga-Luttinger (TL) model [7] as the most 
precise low energy description of ID metals. Here we 
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shall not focus on various very characteristic power laws 
predicted by this model, also regarding spin properties 
[8] coming from the charge sector, but rather focus on 
the question how interactions influence the length Ar. 

As a second striking property the TL-model predicts 
charge-spin separation which, interestingly and contrary 
to statements in the literature [8,9], is not spoiled in semi- 
conducting quantum wires unless spin-orbit coupling is 
not exceedingly strong r\ > 1 to devastate the effective 
mass description. On the other hand, for non-quadratic 
dispersion relations, charge-spin separation is in general 
destroyed with H so . An example are carbon nanotubes 
where v± — vp ± a with a originating in this case from 
the curved surface [10] so that v + ^ v_ spoils charge- 
spin separation. 

A realistic interaction in ID 



V(x) 



1 



1 



K V Vx 2 + w 2 \Jx 2 + w 2 + 4i? 2 



(4) 



is specified by the wire width w and a cut-off at large 
distance R due to screening by the nearest metals, such 
as the gates that help defining the quantum channel. 
Since spin properties involve finite momenta q of the 
Fourier transformed interaction V(q) (and contrary to 
charge properties not the dominant long wave length 
limit V(q — > 0)) the interaction range R plays only a 
minor role if (R/w,k F R) 3> 1. In carbon nanotubes, 
for example, interactions will alter the spin velocity only 
weakly, in marked contrast to the charge velocity [11]. 

How to include the Rashba term in the TL-model ? 
In previous work [8,9] the Fermi velocities v+ and V- 
have simply been put to different values, which for 
e(k) = k 2 /2m* does not describe the leading effect of 
_ff so . Rather both velocities change slightly but obey 
v + = u_ and charge-spin separation. In effective mass 
systems H so acts solely in the spin sector of the corre- 
sponding TL low energy model (of length L), where the 
topological term 



— (v N N 2 + VjJ 2 ) - m*a*v F J a 



(5) 



is most relevant for the following. iV CT and J a denote 
the usual currents of velocities v N/J where with Coulomb 
repulsion (w N ,fj) ^ v F [7]. In strictly spin isotropic sys- 
tems v N — Vj . Since we expect this isotropy to be broken 
only weakly, both of these velocities should be similar in 
magnitude and also similar to the velocity v a of spin wave 
propagation. This latter quantity has been determined 
recently by extensive quantum Monte-Carlo simulations 
[12]. With increasing interaction strength, equivalent to 
a decreasing carrier density, v a /v F was found to decrease, 
even below 0.5 when accounting for parameters of exist- 
ing quantum wires [5] . 

Many quantities of interest can be calculated exactly 
using (5). In particular it can be shown [13] that spins 
polarized in x-direction along the wire precess now over 



a length that acquires a factor v,,/v F compared to the 
length of Eq. (2). With Vj = v a we conclude from 
[12] that this length decreases with increasing interac- 
tion strength. A similar conclusion has been drawn for 
two-dimensional electrons after treating the interaction 
perturbatively [3]. This trend is opposite to what is ex- 
pected for almost linear single electron dispersions (such 
as narrow gap semiconductors) but agrees with experi- 
mental observations [14]. 

For very large and very small k F w asymptotic expres- 
sions can be given for v a which result in the following 
dependencies for the Rashba length: 



m*a* 



:[l-/(2A*t 



if f(2k F w) < — 
Ztt 



at large k F} f(x) 



2- -3/2, 



and 



Af 



tt 2 a-Q 2k F w 
3m* a* w ln2R/w 



if k F <C 1/R 



as is the Bohr radius. This latter estimate follows from 
the conjecture [15] that the system falls into the uni- 
versality class of the Hubbard model at small fillings in 
which case the interaction range R becomes relevant. 

Now we turn to the case of two occupied channels in 
the quantum wire. Then the renormalizations a — > ctj 
and m — > mj due to H so acquire a channel index j; for 
example, in a harmonic confining potential 012 of the up- 
per channel renormalizes by a factor of 3 more than ct\. 

Between the channels acts the same microscopic 
(screened Coulomb type) interaction (4) as within each 
channel. Up to now no theory is available to determine 
spin velocities in multi-channel situations and we would 
like to sketch a possible solution here [?]. We assume 1) 
sufficiently different carrier densities in both channels to 
conserve the TL-phase [16] (otherwise spin gaps could 
spoil the Rashba precession), and 2) charge-spin separa- 
tion. The Fermi momentum in the upper channel may 

be taken as k 2 = ^Jj^kf — 2m 2 wo [17] where cj q is the 

subband spacing and m 2 /i are the effective masses of the 
upper/lower channel after renormalization through H so . 

We expect two spin velocities v a \ and v a i which both 
will be relevant for the Rashba precession, similar to the 
single channel case. They can be obtained from the mi- 
croscopic interaction through generalized susceptibilities 



L d 2 E Q 
tt dNidNj 



and 



, > L d 2 E 
(t»j)« = -- 



7r dJidJj 



(7) 



when suitably generalizing the relation v, 
served in one cha 
(5) in channel j. 



tJv n Vj ob- 
served in one channel. Here, Nj and Jj are the currents 
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How can we obtain v K and Vj for a given microscopic 
interaction ? In principle we know already from the sin- 
gle channel case how notoriously difficult spin velocities 
are to evaluate [12]. Leading corrections at not too small 
carrier densities to the ground state energy Eq/L per 
length can be obtained perturbatively. Fock-exchangc 
terms ~ V(ki ± kj) have to be included adequately for 
spin properties, which goes beyond the RPA approxima- 
tion mostly used for estimating charge velocities [18] since 
the V(q = 0) contribution drops out here. 

The result is 



— V+ V2 



-v+ 

V(2k 2 ) 



V- 



and 



vi + V- —V- 

-V- V 2 + V- 



(8) 



(9) 



Here, Vj are the bare Fermi velocities and V± = (V(k\ — 
fc 2 ) ± V(h + fe 2 ))/27r. Both limits, fc 2 -> (almost 
empty upper channel) and fc 2 — > fci (limit of equiva- 
lent subbands), indicate instabilities in v N and v 3 when 
V(q = 0) > -KVj (long range interactions, note that 
V{q) decreases with increasing q for any realistic electron- 
electron interaction). These are precursors of the Cooper 
or charge density wave instabilities [16], that occur 
in repulsively interacting two-channel systems near the 
threshold for opening the second channel or near equal 
carrier densities in both channels, respectively. We also 
see in (9) that v 3 leaves the Galilei mode (1,1) indepen- 
dent of the interaction. For v\ = u 2 this mode is an 
eigenmode so that perturbation theory pretends Galilei 
invariance of the spin sector, similar as in the single chan- 
nel case [15]. Then perturbative renormalization of the 
spin velocities is solely due to v N . 



> 




kj w 

FIG. 1. Spin velocities in a two channel quantum wire in 
units of the larger of the Fermi velocities vi, versus kiw in 
perturbation theory. 



Resulting spin velocities are shown in Fig. 1. One of 
the velocities, v a i, is decreasing with decreasing carrier 
density 2k\/-K, similar as in the single channel case [12]. 
More strikingly, w ct2 first increases and exceeds the larger 
of the two Fermi velocities v\. Then according to Eq. (2) 
we have evidence for a suppression of the Rashba pre- 
cession as a result of inter-channel coupling by electron- 
electron interaction. 



In conclusion, we have discussed how a hard or soft 
confining potential, and the intra- and inter-channel in- 
teraction affect the Rashba precession along a quantum 
wire. 
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